SQUEEZING IN A 2-D GENERALIZED OSCILLATOR 
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Abstract 



A two-dimensional generalized oscillator with time-dependent parameters is considered 
to study the two-mode squeezing phenomena. Specific choices of the parameters are used 
to determine the dispersion matrix and analytic expressions, in terms of standard hermite 
polynomials, of the wavefunctions and photon distributions. 

o 

1 Introduction 

In the middle of the sixties and beginning of the seventies a set of quantum states of the elec- 
tromagnetic field were observed which have less uncertainty in one quadrature than a coherent 
state [1-3]. These one-mode squeezed states have generated big expectations in optical commu- 



nication systems [4]. In some quantized fields, the interaction hamiltonians occur only between 
pairs of modes and then to understand the main features of the system, one restricts to study 
one and two normal modes . In the last decade two-mode squeezing phenomena have attracted 
attention to study properties of noise and correlations [5-8]. Recently the accidental degeneracy 
of a two-dimensional (2-D) harmonic oscillator with frequency Uq plus an interaction proportional 
to the z-th projection of the angular momentum was studied [9]. This system was called the 
generalized 2-D harmonic oscillator because presents a bigger accidental degeneracy depending on 
the strength A of the angular momentum interaction. This model was generalized [10] to include 
time-dependent parameters, m = mof(t) and A = u \ (t). If we take /(0) = Ao(0) = 1; the 
hamiltonian, for t = 0, represents a charged particle moving in a constant magnetic field . 

The aim of this work is to study two-mode squeezing phenomena with this model because it 
demonstrates the change of dispersions due to variation of the mass and coupling constant during 
the evolution. In the framework of quantum optics the hamiltonian is built by: the operator 
(1// + f)a) • a, that causes a time-dependent exchange of kinetic and potential energies within 
each mode; the interaction 1/2(1// — f)(cP ■ eft + a ■ a), which describes a degenerate two-photon 
interaction; and the potential i\ (t)(al ), that is a mode mixing operator. 

The solution of the corresponding time dependent Schroedinger equation is obtained through 
the theory of integrals of motion [11]. By means of Noether's theorem, using a special variation 



[10], we construct the linear time dependent integrals of the motion. The resulting quantum 
invariants are given in terms of the positions and momenta operators [10,11] by 



P(t) = X,p + X 2 q, Q(t) = X 3 p + X 4 q , (1) 

with the initial conditions P(0) — p and Q(0) = q, so that the 2x2 matrices previously introduced 
satisfy Ai(0) = A 4 (0) = I 2 and A 3 (0) = A 2 (0) = 0. The operators A(t) = l/y/2 [Q(t)/l + il/hP(t)] 
and its hermitean conjugate, can be constructed with the matrices 

Ap = yA 3 + — Ai, \ = yA 4 + — A 2 , (2) 



with / = \/ rap^p defining the oscillator length. These integrals of motion also are given in terms 



of the creation and annihilation photon operators 

A(t) = M x a + M 2 a\ A\t) = M 3 a + M 4 a ] ■ (3) 

With the initial conditions A(0) = a and A^(0) = eft, the matrices defined in (3) comply with 
Mi(0) = M 4 (0) = I 2 and M 3 (0) = M 2 (0) = 0. The A fe 's, M fc 's, A p and X q are entries of symplectic 
matrices in four dimensions because the invariants (1) and (3) satisfy the commutation relations 
of Heisenberg-Weyl algebras. 

In the present work, we study the behavior of the model for Ao(t) an arbitrary function of time 
and considering two kinds of varying masses, i.e., two choices for the function f(t), namely: 

/(f)=exp(7t) ; (4) 

( 1 , * < 

f{t) = I cosh 2 tt t , < t < T . (5) 

I {Q (t - T) sinh Q T + cosh fi T} 2 , T <t 

For these two cases the A^ matrices take the general form 

Afe = ^ R = ^(-sl1 co n s^) ; fc = 1 ' 2 ' 3 ' 4 ' 

where the definition 9 = J f uj \ (T)dT was used. The analytic expressions for the //^'s functions 
are given in Ref. [10]. In the next sections we determine the coherent and Fock-like states, the 
photon distributions and the dispersion matrices in terms of these /i^'s. 



2 Squeezed Coherent and Fock States 

The coherent-like states are obtained by solving the differential equation A(t)Q (q,t) = with 
A(t) given in Eq. (3). This solution yields the vacuum state of the physical system, and its phase 
is chosen to guarantee that satisfies the time dependent Schroedinger equation. The expression 
for the ground state wavefunction is 



To get the last expression the relation (2) was used and the functions n P = -j= (j-/J>i + 7^3) and 
= 73 (f^ 2 T'" 4 ) were defined. To obtain the general expression for the eigenstates in the 
coordinate representation one needs to apply the unitary operator D(a) = exp{c? • — a* ■ A}, 
which is an invariant, to the vacuum wavefunction (7), i.e., 

$ a (q,t) = exp|-^ + ^a-a + ^-gRa|$ (g,t) . (8) 

These are expressed in terms of mult i- dimensional Hermite polynomials [12] through the relation 

Substituting the last expression into (8) and using the form of the coherent-like states in the 
Fock-like representation, we get the Fock-like eigenstates in the coordinate representation: 



(gjnma) = $ (q,t) H„\,„; > ['J^W) ■ (10) 

These multi-dimensional Hermite polynomials are rewritten as a product of two standard one- 
dimensional Hermite polynomials [12] as follows: 



H - 1 (-^) - (-£) H " i vm^\ |cos " 91 + 81,19 *0 

x H n2 ( ^_ * [- sin 6 q x + cos 6 g 2 ]) , (H) 

where we use the explicit expression of matrix R. These Fock (10) and coherent (8) -like states 
represent squeezed and correlated eigenstates of the system as it will be shown further. 



3 Propagator 

The propagator in the coherent state representation is given by the matrix elements of the evolution 
operator U(t), which will be obtained by means of the theory of time dependent integrals of 
motion [11]. If I(t) is an integral of motion then satisfies I(t)U(t) = U(t)I(0). Taking its matrix 
elements with respect to the coherent states, we get a linear system of differential equations, which 
can be solved. Thus the propagator takes the form 

G(a*,7,t) = ^ H fJJ Mi ^ e X p(-^*Mr 1 M 2 ^ + ^Mr 1 7+ \lM*M?j) . (12) 
For the cases (4) and (5), the following relations are satisfied 

>JdetM 1 = -±= (l^ - = g 1} Mf 1 M 2 = -±- (l N + j ^ I 2 = g 2 I 2 , (13) 



M^ 1 — ^-R , M 3 M 1 - 1 = ^-(//i; + j^I 2 = g 4 h ■ (14) 

Substituting these relations into the Eq.(12) we get the propagator, which through Eq. (9) can 
be expressed in terms of multi-dimensional Hermite polynomials. If we compare with the power 
series expansion of the propagator we get the probability amplitude for having rt\ and n 2 photons 
in the coherent-like state I7, t), i.e., 

( ™'* t>= i^r p ( J f ■ (15) 

By means of the Eq.(12) this amplitude can be rewritten in terms of standard Hermite poly- 
nomials [12]. The squared absolute value of this amplitude yields the photon distribution function 
of the system, W nin2 (^,t) = |(nin 2 |7, t)\ 2 . This will let us calculate, at least formally, the mean, 
(N k ), and the mean squared fluctuation of the number of photons, (AiV fc ) 2 , in direction k, which 
are present in the coherent state (7, t). The expectation values of N k and N% are evaluated directly 
using the expressions of the creation and annihilation photon operators in terms of the integrals 
of the motion (3), and the commutation properties for these invariants. For the vacuum state one 
has 

(N k ) = 7 (jui -/U 4 ) 2 + (m uJoH3 H A*2 ) \, (16) 

( N k) = ^{(^1- + (m uJofi 3 + —^—fi 2 ) 2 } 
8 I m uJo J 

| (A*i + At4) 2 + (™o^o^ 3 —fi2) 2 \ + 777 ((/"I - Ak) 2 + {moUJo^ H —^ 2 ) 2 \{17) 

L m uj J 16 I m uj J 

With these expressions, we evaluate the ratio of the mean squared fluctuation (AN k ) 2 and the 
mean number of photons (Nk), which determines the nature of the distribution function of the 
system: 

(AiV fc ) 2 If, v 2 / 1 \ 2 1 

(N k ) 2 I V m uj ) J 

For the cases (4) and (5) the ratio is greater than one when t > 0, which implies that we have 
a super- Poissonian photon distribution function. For t — 0, there is a discontinuity in the ratio, 
which is obtained by comparing the following limiting procedures: making t — > and then a — > 0, 
and conversely. 



4 Dispersion Matrices 

The dispersion matrix can be written in terms of 2 x 2 matrices characterizing the dispersions 
in the positions and momenta operators and the correlation between them. Besides for the cases 
under study, due to (6), they take the form 

= l hm ° Uo (j^o)~^ 1 + ^) 12 ' (19) 



°"?pCO = -%{—— A*iA*2 + mouo^/j,^ I 2 , (20) 

<^(*) = ij (/i? + (m ^ ) V3) I2 • (21) 

The corresponding correlation matrices for the creation and annihilation operators are obtained 
immediately from the last expressions; they are given by 

al = \{£ - & + (mo^o) 2 ^ " j^-y - 2. (^/^ + } , (22) 

2 1/2,2,/ \2 2 , ^2 m 0^o\ / OQ x 

a a t a = i ^ + /x 4 + (mo^o) /i 3 + (mQL , o) 2 j > ( 23 ) 

^otat = \ [lA - ^4 + {m Q U ) 2 l4 - + 2i (^^^ l/i2 + m 0^0/i3^) J • (24) 



Fig. 1. Dispersion and correlation matrices behavior in positions and momenta 
space for the studied cases in this paper: (a) corresponds to Eq.(4), and (b), to Eq.(5). 



The behavior of the dispersion matrices is illustrated in Fig. 1. For the case (4), we choose the 
parameters 7 = 0.1 and mo = ujq = 1. It is seen that there is squeezing for the coordinates and 
stretching for the momenta. Also one notes that a pq is a negative function and therefore there 
are one-mode correlations between the coordinates and the momenta. If we reverse the sign of 7, 
the roles between the dispersion for coordinates and momenta are interchanged, and a pq becomes 
positive. In the case (5), we use the parameters VLq = 0.15, T = 10, and mo = ujq = 1. In spite 
of the mass is different that in the previous example, the general trends are similar. For example, 
the a pp is an increasing function of time starting from its minimum value at t < 0, and there is 
squeezing for the a qq . The main difference appears in the correlation a pq : in this case, it can be 
positive for large times, while in the previous one is negative or zero for any time. 
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